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Abstract. Given a real vector space V of finite dimension, together with a 
particular homogeneous field of bivectors that we call a field of projective forces, 
we define a law of dynamics such that the position of the particle is a ray i.e. 
a half-line drawn from the origin of V. The impulsion is a bivector whose 
support is a 2-plane containing the ray. Throwing the particle with a given 
initial impulsion defines a projective trajectory. It is a curve in the space of 
rays <S(V), together with an impulsion attached to each ray. In the simplest 
example where the force is identically zero, the curve is a straight line and the 
impulsion a constant bivector in /\ 2 V. A striking feature of projective dynamics 
appears: the trajectories are not parameterized. 

Among the projective force fields corresponding to a central force, the one defin- 
ing the Kepler problem is simpler than those corresponding to other homo- 
geneities. Here the thrown ray describes a quadratic cone whose section by 
a hyperplane corresponds to a Keplerian conic. An original point of view on 
the hidden symmetries of the Kepler problem emerges, and clarifies some re- 
marks due to Halphen and Appell. We also get the unexpected conclusion that 
there exists a notion of divergence-free field of projective forces if and only if 
dim V = 4. No metric is involved in the axioms of projective dynamics. 

1. Introduction 

The main object of the present study is the class of differential systems defining, 
on an open set U of an affine space A, the motion of a particle subjected to a 
field of forces. We write such a system: 

q = f(q), qeUcA, f : U - A. (1.1) 

The force or acceleration f(q) lives in the vector space A associated to A. An 
elementary computation, apparently due to Appell, indicates that it makes sense 
to "projectivize" such systems, considering the n-dimcnsional affine space A as 
an affine hyperplane of a n + 1-dimensional vector space V. Then A plays the 
role of an affine chart for V(V), the real projective space associated to V, or 
better for S(V), the double covering of V(V), whose points are the rays from the 
origin O of V. The topological aspects of this projectivization are interesting, 
but quite trivial: one can sometimes extend the domain U of the motion; for 
example one "closes" the hyperbolic orbits of the Kepler problem, making them 
ellipses. The local aspects are more surprising: Appell's computation indicates 
that something as (1.1) is already defined when the particle lives in a space which 
is not affine and whose tangent bundle is not endowed with a linear connection. 



1.1. Appell's computation. It may be presented as follows. As A C V, the 
"point" q E A is now a "vector". Let /i E V* be the linear form such that 
(h, q) = 1 is the equation of A. We choose a non-zero hi E V™ and call Ai 
the affine hyperplane of equation {h\,q) = 1. To a g £ 4 with (hi,q) > we 
associate gi = (h\,q)~ l q E Ai on the same ray. We compute 

. = (hi,q)q - (hi,q)q 

Because of the denominator, the second derivative q\ is quite complicated. But 
we change the time parameter. We define the derivative r' of a quantity r with 
respect to the new time by 

r' = (hi,q) 2 r. (1.2) 

Then 

q'^ihuqfUh^q^-ihuqjq). (1.3) 

Substituting q — f(q) and q = (h, qi)^ 1 qi, we observe that the right hand side, 
which is a vector tangent to A\ at q%, only depends on the position q\. Thus 
a system such as (1.1), i.e. defined by a field of forces depending only of the 
position, remains of the same type after a "change of projection" , provided the 
time parameterization is changed according to the rule (1.2). 



h = 1 




'O 

1.2. Figure 

1.3. Halphen's computation. It particularizes the previous computation and 
introduces the "classical gravitation" aspect. We consider a center of force 
c £ A. We consider the particular case of (1.1) 

q=-\\q-cf(q-c), qeA\{c}, € H. (1.4) 

We endowed the affine space A with a Euclidean structure ||.|| 2 , making it a 
n-dimcnsional Euclidean space. We apply the above transformation 

q'{ = -\\q-cf(h 1 ,q) 2 ((h 1 , q )(q-c)-(h 1 ,q-c)q) 
= -\\q-cf{h 1 ,q) 2 (-(h 1 ,q)c+(h 1 ,c)q) 



or 



3i 



|g-c||^(^,g) 3 (^,c)( t^-t-t^-t 
\{hi,q) (fti,c) 



We shall express the right hand side in term of q\ — (h\,q) 1 q E A\ and 
ci = (hi,c)~ 1 c E Ai. For this we extend the Euclidean quadratic form ||.|| 2 
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from A to V in a special way. We decompose any q G V in horizontal plus 
vertical components: q = qA + Ac, where qA € A and A € R. Here we think 
of A as the horizontal vector hyperspace in V, of equation (h, q) = 0. We set 
Ikll 2 = Ik^H 2 ; in particular ||c|| 2 = 0, and \\q - c|| 2 = ||q|| 2 . We resume the 
computation: 

q'l = -\\qf{h 1 ,qf{h 1 ,c){q 1 -c 1 ) = - \f(h u qf +fi (h u c)( qi - Cl ). 

As (h,qi) = we get 

«i = -Iki-crf^g!)- 3 "^^^)^ - Cl ). (1.5) 

The transformed system (1.5) has the same form as the original system (1.4) if 
/3 = —3, i.e. if (1.4) defines the Kepler problem and corresponds to the Newto- 
nian attraction. Actually we must slightly extend the form (1.4) to accept the 
"mass" factor (hi,c) in (1.5). Moreover, the affine hyperplane A\ is endowed 
with the restriction of the degenerate quadratic form ||.|| 2 to it. If A\ is "verti- 
cal" , i.e. if it contains the vertical vector c, this restriction is degenerate; again 
we shall slightly extend the class of equations (1.4) to accept this kind of ||.|| 2 . 

It is important to have clearly in mind that if (3 ^ — 3 the transformed system 
(1.5) is more complicated than (1.4). 

1.4. Divergence of the field of forces. The next statement could be obtained 
by a long but straightforward computation. We will prove it simply in §5. If 
dim A = 3 and div/ = in (1.1), then the field of forces in the right hand side 
of (1.3), expressed in the new variable q\, is also divergence- free. 

This statement is readily checked in the example of the Kepler problem. It 
is well-known that in dimension 3 the force function in this problem satisfies 
the Laplace equation, and that consequently its gradient, the field of forces, is 
divergence- free. 

The Kepler problem has another striking property: its bounded orbits are pe- 
riodic. Newton emphasized this property. For him one tests the inverse square 
law verifying that the aphelia of Saturn and the comets are at rest. This law is 
(1.4) with (3 = —3; the force field is divergence-free. Is this coexistence between 
the dynamical fact and the null divergence a mere coincidence? We do not know 
if this question has already been discussed. Maybe it seems too vague, especially 
if we stick to the unique Kepler problem and do not discuss other examples. In 
the 19th century, examples with both the dynamical and the divergence-free 
properties were discovered. They are discussed in [Alb], [A1S] and [BoM]. We 
will present the most general example we know in §6. 

Does projective dynamics explain this coexistence? Choosing the projective 
point of view we "forget" the affine and metric structures of the space. But we 
keep a structure which is sufficient to express the divergence-free property, the 
central force property and the degeneracy of the dynamics. So, if some expla- 
nation is to be found, it is reasonable to look for it inside projective dynamics. 
Actually when we think that there exists a relation between these properties, 
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we immediately raise two objections. The first one remains for us mysterious: 
the dynamics is two-dimensional, while the divergence-free hypothesis is three- 
dimensional. 

A second objection is rather a question: does a field of forces that is divergence- 
free in dimension n define a special dynamics? Or is there something special only 
if n = 3? A deep property of the case n = 4 was discovered by Jacobi (sec [A1C] 
pp. 161 and 169, [Mon]). Concerning the central force field which is divergence- 
free in dimension n = 2, we do not know anything remarkable, except that 
Johann Albrecht Euler, son of Leonhard Eulcr, obtained an expression for the 
total duration of a rectilinear free fall, with a zero initial velocity. As mentioned 
to me by James E. Howard, this result is an elementary exercise requiring the 
computation of J °° e~ v / 2 dv. 

Projective dynamics offers a surprising answer to this second objection: in this 
"poorly structured" context, the concept of a divergence-free field of forces only 
exists in dimension 3. 

2. s-tangent vectors and s-scalars 

As above, V is a real vector space of dimension n+ 1. A ray is an open half-line 
drawn from the origin in V. In projective dynamics, the position of a particle 
is a ray. 

We call S(V) the manifold of the rays of V. To describe in words and symbols 
a relation between V and S(V), we use such expressions as: a point Q G S(V) 
corresponds to a ray Q C V. Actually there are two presentations of the pro- 
jective facts. One is n + 1-dimensional and elementary; it corresponds to the 
classical homogeneous coordinates. The other one is n-dimensional, which is 
more natural; but it leads to abstract constructions. 

2.1. Definition. Let s G H. A s-tangent vector w to S(V) at Q G S(V) 
corresponds to a class of homogeneous vector fields of degree s along the ray 
Q C V. We define a class w as follows. Let w : Q — > V be a homogeneous 
vector field in w. Then wi : Q — > V is in the same class if and only if the vector 
field wi — w is tangent to the ray Q. 

Given an open set U C S(V) we denote by TqU the n-dimensional vector space 
of s-tangent vectors at Q G hi. The description of a tangent vector to P(V) 
or to S(V) as a 1-tangent vector is standard. It is the identity TqU = TqU. 
For the values s ^ 1 of the homogeneity, a non-zero s-tangent vector is not a 
tangent vector, but nevertheless it points a direction tangent to U. 

2.2. Equivalent definition. Let s e H. A s-tangent vector w to S(V) at 
Q G S(V) corresponds to a w : Q — > f\ 2 V satisfying q/\u(q) = for any q G Q, 
and w(Xq) = \ s+1 w(q) for any A > 0. 

We defined w first as a class of homogeneous vector fields along the ray Q C 
V, then as a homogeneous bivector field along Q. It is easy to relate both 
definitions. If wq '■ Q — > V is in the class w, we associate to it the bivector field 
qAw (q), which is homogeneous of degree s+1 and satisfies qAqAw (q) = 0. If we 
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take another vector field w\ in the same class, then qAwg = qAwi. Conversely 
we consider w : Q — > f\ 2 V and we choose any q e Q. As q A w{q) = there 
exists a v e V such that w(g) = q A i>. We set w n (q) = v and extend w by 
homogeneity. 

The notion of a s-tangent vector is extremely useful in projective dynamics. 
This is why we introduced it first. Nevertheless the following definition of a 
generalized "scalar" quantity should logically come first. 

2.3. Definition. Let set A s-scalar p at Q £ S(V) corresponds to a function 
p : Q — > H such that p(Xq) = X s p(q) for any q £ Q and any A > 0. 

A 0-scalar is simply a real number. For a U C 5(F) we denote by 9 S U the 
line bundle (i.e. the vector bundle with one-dimensional fiber) of s-scalars. The 
notation 0(s) for the same object is widely used in the case of a complex 
projective space, s being an integer. The following formulas are quite standard: 

e s u = ^u)® 8 , e r u ® e s u = e r+s u, e s u ® t t u = T r+S u. 

Here the tensor products of fields defined on S(V) correspond to mere products 
of the corresponding fields defined on V. 

3. The data used in projective dynamics 

3.1. Definition. We call a (— l)-tangent vector a projective impulsion and a 
(— 3)-tangent vector a projective force. 

A field F of projective forces onUc S(V) is a section of the vector bundle 
T~ 3 U. With a standard notation this reads F G F(T~ 3 U). To describe the 
"corresponding" object, we denote byWeF the union of the rays corresponding 
to the points of U C S(V). The field F corresponds to a map 
positively homogeneous of degree —2 and such that q A F(q) = at any q <EU. 

A parameterized path is a (smooth or real analytic) map (p : I — > U where 
I C H is an open interval. To define an oriented (unparameterized) path, we 
simply weaken the structure of the source space I. 

3.2. Definition. An oriented path is a map c : 1 — > where 1 is a 1-dimensional 
oriented manifold diffcomorphic (or analytically diffeomorphic) to H. 

While working with autonomous systems as (1.1) it is natural to consider that 
a solution is more than an oriented path, but less than a parameterized path. 
We mean we do not consider that the scalar value of the time t at a point of the 
path is a relevant information; but the time At to go from a point to another 
is well-defined. We call such solution a usual trajectory. 

It is useful for our purpose to describe a usual trajectory in the following way: it 
is an oriented path together with a field of velocity vectors. The velocity vector 
at a point must be tangent to the path at this point. To take into account 
the possibility of singularities and multiple points, we should endow the source 
manifold X, rather than the image of c : X —> U, with a field of non-zero, 
positively oriented vectors. Such data induces a map X — > TU. 
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A projective trajectory is an oriented path together with a field of tangent 
projective impulsions. This is the main idea. However the possibility of multiple 
points and singularities somewhat complicates the precise definition. 




3.3. Figure. A usual trajectory and a projective trajectory. An attempt of 
visualization. 

3.4. Definition. Let c : X — > U be an oriented path. Let OX be the line bundle 
above X, pull-back by c of 9~ 2 U. A structure of projective trajectory on c is a 
positive section of the oriented line bundle TX ® OX. A projective trajectory is 
an oriented path together with a structure of projective trajectory. 

3.5. The exponent —2 in the definition of the bundle OX is the difference 
between the homogeneity of a projective impulsion and the homogeneity of a 
tangent vector. To put this abstract definition into practice, let us show that 
a projective trajectory (c, a), where c : X — > U and a G r(TX<g> OX), induces 
a map to the projective impulsions c a : X — > T~ X U. At x G X we choose any 
non-zero p G T*X, which determines a unique v G T X X such that (p,v) = 1. 
The contracted product p ■ a\ x G O x X defines a (— 2)-scalar A G S~ 2 U : where 
y — c(x). Let us call c*v G TylA the push-forward of v. The resulting projective 
impulsion c a {x) = A ® c*u G Ty X U does not depend on the choice of p: if we 
multiply p by a non-zero real number a, we multiply A by a and divide v by a. 




3.6. Figure. A projective trajectory and a projective impulsion. 

As always the definition of the corresponding object is more elementary. To an 
oriented path c : X — > hi corresponds a two-dimensional "half-ruled" submani- 
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fold of X x IA. By half-ruled we mean that it contains half-lines instead of lines. 
The projection of this object on the factor IA is a two-dimensional semi-cone, 
possibly with multiple points and singularities. A structure of projective tra- 
jectory is a field of bivectors tangent to this semi-cone, possibly many-valued. 
The "tangency condition" must take into account the possible singularities of c. 

3.7. In the definition below and everywhere in this text, X is an oriented manifold 
diffcomorphic to HI. We denote by / an open interval, and r : I — > X, t r(t) 
a global chart of X, which respects the orientation. 

3.8. Equivalent definition. A projective trajectory corresponds to an oriented 
path c : X — > U together with a map 7r : X — > /\ 2 V satisfying the following 
"tangency condition" : there exist a global chart r : I — ► T and a map q : I — > V 
satisfying q(t) e Q{t), where Q(t) — c(r(t)), and n(r(t)) = q(t) A q(t). 

3.9. To show that both definitions are equivalent, we start from a projective 
trajectory (c, a) in the sense of Definition 3.4. We construct c a as in 3.5 and set 
7T = c a . We shall check that ir satisfies at any x £ X the "tangency condition". 
Let t ; / — > X be a chart 3.7. We set r(t) = x, Q(t) — c(x), and construct with 
p = dt a A e Sq 2 U as we did in 3.5. Then c*v corresponds to a civ : Q — > /\ 2 
such that c*t;(g(t)) = A q(t) as soon as G We select the unique 
q(t) in the ray Q(t) such that X(q(t)) = 1. The tensorial product relation in 3.5 
means that 7r(x) = c*v(q(t)), which gives the result. The same computations 
allow to construct the section a if we start with Definition 3.8 of a projective 
trajectory. 

3.10. Why did we choose s = — 1 for the projective impulsion? We will see 
that it is forced by our introductory remarks. But we have already a strong 
argument to present. Among the projective trajectories we can distinguish 
the linear ones, which "draw" part of a projective line in S(V). Among these 
rectilinear motions we should be able to distinguish the "uniform motions". 
The statement "constant projective impulsion" sounds good, and it makes sense 
with our choice s = — 1. In this case, by the "bivector definition" of a s-tangent 
vector, the bivector is constant along the rays. To define the uniformity of the 
motion, we require that the bivector is constant on the whole vectorial plane of 
V that corresponds to the rectilinear trajectory in S(V). 

4. Prom screen dynamics to projective dynamics 

In this section, the statement that occupies us is: in a given field of projective 
forces, a "position ray" and a projective impulsion at this ray uniquely determine 
a projective trajectory. To establish it, we study the process that converts our 
quite abstract "projective dynamics" into a more familiar "screen dynamics" , 
whose simplest case is (1.1). The converted statement is the usual statement 
that an initial position together with a velocity uniquely determine a trajectory. 
We convert again this trajectory into a projective trajectory. Finally we prove 
that the resulting projective trajectory is independent of the screen used in its 
construction (except if the screen is too small, i.e. if it covers only part of the 
projective trajectory). 
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4.1. Definition. Let U C S(V) be an open set and U C V the corresponding 
open semi-cone. A screen function h for U is a positive /i : W — > IR, positively 
homogeneous of degree one 1 . The associated screen Uh for U is the hypersurface 
of U with equation h = 1. 

The simplest screen function is a linear form on V. It gives a "flat" screen for 
an open "hemispheric" U C S{V). In §1 we used a pair (A, Ai) of flat screens. 

4.2. Restriction of s-tangent vectors to a screen. Consider the screen h = 1 for 
W C <S(V). A s-tangent vector v at Q £ U defines a unique vector Vf l (q) e V, 
tangent to the screen at the unique q <E Q with h(q) = 1, by the formula: 

v h (q) = dh(q)}v(q), (4.1) 

where J is the interior or contracted product, and v : Q — > f\ 2 V corresponds to 
v. The formula (dh(q), Vh(q)} = shows that Vh is tangent to the screen. 

4.3. Restriction of projective trajectories. The points of U C 5(V0 are in one-to- 
one correspondence with the points of the screen Uh C V with equation h = 1 . 
So an oriented path c : T ^ U restricts into an oriented path Ch : X — > . If c 
possesses a structure of projective trajectory, i.e. a map t) : I — > /\ 2 V with the 
tangency condition 3.8, is endowed with the field of tangent vectors it is 
a "usual trajectory" . 

4.4. Law of areas. If it has no multiple points, a projective trajectory is simply 
a two-dimensional semi-cone endowed with a tangent bivector field v, which 
is constant along the rays. The field v defines a way to measure the oriented 
area of any domain delimited on the semi-cone. One takes as unit the area 
of a tangent parallelogram spanned by two tangent vectors a and fj such that 
v = a A f3. 

The usual trajectory Ch on a screen Uh, restriction of the projective trajectory, 
is described according to the "law of areas" : the moving point sweeps out from 
the origin of V equal areas in equal times, more precisely it sweeps out a unit 
of area in two units of time. Indeed, in a time dt the area swept is half of the 
area of the parallelogram spanned by q and Vhdt. And the area of the tangent 
parallelogram spanned by q and Vh is one. 

4.5. The example of the uniform motion. In 3.10 we stated that a uniform 
projective trajectory corresponds to a vectorial plane in V with a non-zero con- 
stant area element v e f\ 2 V, whose support is the plane. From the law of areas 
the reader will deduce that the restriction of this projective trajectory to a flat 
screen is a usual uniform motion, and its restriction to a screen corresponding 
to a screen function h = \\q\\ is a uniform (geodesic) motion on the unit sphere 
H^ll = 1. In this last case we need a positive definite quadratic form q m ||g|| 2 . 
Of course the case of an indefinite quadratic form is also interesting. 

4.6. A field of forces on a screen. Given a field of projective forces onWc S(V) 
and a screen function h : U — > H, a usual field of forces may be defined by 

1 Thus it corresponds to a positive section of the 1-scalar bundle 8 1 U. There is no special 
reason to choose one as the degree of homogeneity. We simply need a convention. 
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restriction to the screen Uh C V with equation h = 1. It is tangent to the 
screen. What kind of dynamics is defined by such a field of forces? 

Let fh be any vector field tangent to Uh- As a force Geld, it defines a dynamics 
through the system 

q = h(q) + Xq, (4.2) 

where the real value of X is forced by the constraint h(q) = 1. The term Xq is 
a "reaction" . "Normal" reactions are more traditional, but in this framework 
we have to take a "radial" reaction. The determination of A is a valuable 
exercise. As h(q) = 1, (dh(q),q) = and (dh(q),q) + (d 2 'h(q) , q ® q) = 0, 
where d 2 h : V —* V* ® V* is the Hessian quadratic form. By Euler's relation 
{dh{q),q) = h(q) = 1. This gives A = -(d 2 h(q),q® q). 

It is natural to consider that vector fields are "simpler" than bivector fields. 
However, everywhere in this theory bivectors appear as a simplifying tool. This 
happens even while working with screen dynamics. Equation (4.2) above be- 
comes 

7T = /(<?), with ir = qAq, /(«) = g A //,(?). (4.3) 

From this equation we can deduce (4.2): contracting dh at the left we get 
dh\ (q A q) = (dh, q)q — (dh, q)q = dh\f — fh- As (dh, q) = 1, this equation is of 
type (4.2). 

4.7. Two screens. Suppose a field of projective forces is given by a / : U — > A 2 ^ 
positively homogeneous of degree —2. Suppose a screen h(q) = 1 is given, and 
the force is restricted to the screen. Suppose a trajectory q(t) on the screen is 
found, solution of the screen motion equation dir/dt = f(q), with tt = qA(dq/dt). 
Then a projective trajectory is defined extending n(t) on the ray of q(t) with 
the homogeneity zero. 

Suppose qi(t) — X(t)q(t), X(t) e H, is permanently on another screen. If t is a 
new parameter such that dr/dt = [X(t)] 2 , then ni(t) = qi(t) A (dqi/dr) = n(t). 
We have d^/dr = (d% / dt)(dt / dr) = X- 2 f(q) = f (Xq) = f( qi ). Then qi {r) is a 
trajectory on the other screen, solution of the screen dynamics induced by the 
restriction of / to this screen. 

Thus starting with a ray and a projective impulsion, we construct a projective 
trajectory using a screen; this trajectory is unchanged if we use another screen. 

4.8. A more intrinsic equation. Let t i— > (Q(t),n(t)) be a projective trajectory 
on U C S(V), with Q(t) e U and n(t) E f\ 2 V . There exists a vector field 
v along the two-dimensional semi-cone corresponding to t \— ► Q(t), such that 
n(t) = q A v at any t and any q G Q(t) C V. We should think v as the velocity 
q; it is not unique but determined at q up to the addition of Xq, with A £ R. If 
a field of forces / is given as a field of bivectors on U, and if we think 7r as a field 
of bivectors along the semi-cone, the equation of dynamics may be written: 

CU - f(q). (4.4) 

Here d v means the standard derivation in the direction of the field v. Note that 
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this equation coincides with (4.3), that d v+ \ q ir = d v ir + Xd q TT and that d q n = 
because 7r is homogeneous of degree 0. 

4.9. Remark on a notation. When we add v to Xq we implicitly consider that 
they are object of the same type, i.e. vector fields along the semi-cone. We will 
consider q as a notation for the "Eulcr field" , i.e. the vector field on V whose 
value at q is q. Here this field is implicitly restricted to the semi-cone. 

5. Divergence-free fields of projective forces 

To study the divergence operator acting on vector fields, we first need to intro- 
duce the generalized differential forms on U C S(V). We build up an algebra of 
tensors on the s-tangent vectors introduced in §2. We do not fix a terminology 
and try to work with a minimal system of notation. We adopt a simplified 
notation for T°U and forget the notation T S U, using instead: 

XU = T°U, xu <g> e s u = T S U. 

The generalized tensors fields are the sections of (g)* X*U ® J XU ® 9 S U. The 
"usual" tensors fields, i.e. those constructed considering that U has only the 
local structure of a differentiable manifold, are sections of ^ T*U ® ® J TU. 
But TU = T X U = XU ® 9 X U and T*U = X*U <g> B^U; expanding this we get 
(g)* T*U ® <g> j TU = <g) f X*U <g> <S> j XU <g> 9 j - l U. 

A u) e T(/\ l X*U ® 9 S U) corresponds to a lu : U — > f\ V*, positively homo- 
geneous of degree s, and such that q\uj(q) = for any q E U. The following 
identities hold at the point q G U: 

q\duo = C q u = (i + s)cb . (5.1) 

The first is Cartan formula. Let d q be the standard derivation, in the direction 
of the Euler field q, of functions defined on V. We get, using coordinates or con- 
siderations on the derivations, £ q ui — d q uo = id). But by Euler's characterization 
of homogeneity d q ib — sub, which gives the second formula. 

5.1. Proposition. Let oj £ r(/V X*U <S> 9 S U) and let u) : U -» /\ l V* be the 
corresponding i-form. Suppose there exist s' s H and p G r(/\ 4+1 X*U (£> 9 s U) 
such that duj = p. Then i + s = 0: w is a "usual" i-form on U and consequently 
possesses an exterior derivative duj. Finally duj = p and s' = s — 1 = — i — 1. 

Proof. Assuming duo = p, we get q\du = and i + s = by (5.1). Finally 
duj satisfies dili — dw, because " is the pull-back operation by the canonical 
projection U — > ZY, which commutes with d. 

Divergence-free fields of forces. Let h : U — » 1R be a screen function, and 
//j : W/j — > V be a vector field tangent to the screen % with equation h = 1 . If 
we choose a unit of volume on V", an area form is canonically defined on Uh and 
the divergence of fh is well-defined. 

Let n G A" +1 V* be the non-zero volume form defining the unit of volume in 
V. The n-dimensional area form at q G Uh is simply the restriction to Uh of the 
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form gj/i. One can also think of it more geometrically. The area of an open set 
in Uh is n + 1 times the volume of the cone with base the open set and vertex 
at the origin of V. 

The divergence of fh is the function on Uh defined by the classical identity 
(divf h )q}(i = d{f h \{q\n)). We set / = qAf h ; then fhl(ql^) = -f\H- We extend 
/ by positive homogeneity of degree —2 in a projective force field / : U — > /\ 2 V. 
The vector field fh is divergence-free if and only if the restriction to Uh of the 
n-form vanishes identically. 

This last condition is not invariant under a change of screen. However, if d(f\u) 
vanishes identically on a screen Uh, and if moreover q\d{f\u) = 0, then it 
vanishes identically on U. It is only in this case that we can say that the 
projective force field f is divergence- free, i.e. div fh = for any screen function 
h. We can apply Proposition 5.1 to the closed form Co = f\[i, which corresponds 
to a section lo of A™" 1 X*U <g> 9~ 2 U. This gives the following statement. 

5.2. Proposition. Divergence-free projective force fields exist only if n = 3, i.e. 
if dim V = 4. In this dimension, if we choose a unit of volume on V, the space 
of projective force fields onWC S(V) is canonically identified with the space of 
2-forms on U. A divergence-free projective force field is sent on a closed 2-form. 
In this dimension again, if a projective force field restricts to a given screen Uh 
into a divergence-free vector field, it is a divergence-free projective force field, 
and consequently it restricts to any screen into a divergence-free vector field. 

6. Generalized Kepler problem 

6.1. Central forces. Let V be a vector space of dimension n + 1, h € V* a 
non-zero linear form, and A C V the affine hyperplane with equation h = 1. 
Let c£ Abe the "center of force" . Wc define the motion of q G A \ {c} in a 
central force field by the system 

q = -^(q)(q-c^ V : A \ {c} - JR. (6.1) 

Let U C S(V) be the set of rays intersecting A \ {c}, and U the correspond- 
ing semi-cone. We extend by positive homogeneity the force field to obtain a 
projective force field / : U — > /\ 2 V, whose expression is 

/(9) = *(9)?Ac, where ttfa) = (h, g )" 3 v( ) (6.2) 

is positively homogeneous of degree —3 and satisfies ^(q) = tp(q) on A. 

6.2. Let us suppose that V has the homogeneity property of the coefficient 
||<7 — cj | 3 in the Kepler problem; namely 

ip{Q) = <f>(q-c), where tj>{\v) = \~ z cf>(v), A > 0, (h,v)=0. (6.3) 

Then *(<?) = <j>(q - (h,q)c) = <f>(h\ (c A q)). Wc have *(g + 7c) = for 
any (q,j) G U x H such that q + 7c G U. The coefficient is invariant by 
any translation with direction c. Let [c] C V be the line generated by c. The 
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coefficient ^ and consequently the projective force field are extended canonically 
toV\ [c] by this property. 

This canonical extension exists for the Kepler problem (see §7; observe that \P 
is an even function of q.) One can also extend force fields defined by several 
fixed point masses on A, or more generally by any fixed repartition of mass. We 
will only study the situation above, i.e. the motion in a projective force field 

f(q) = V(q)q A C, (6.4) 

where for any q G V \ [c], 7 € H and A > 0: 

*:^\[c]-f]R, *(g + 7 c) = *(«), *(Ag) = X^iq). 

6.3. Lemma. Let G be the group of linear transformations of V that fix the 
vector c and induce the identity on the quotient space V/[c\. An element of G 
sends any projective trajectory solution of Problem (6.4) onto another solution. 

Proof. The projective force field / of (6.4) is trivially invariant by such a 
transformation. 

6.4. Lemma. The bivector "constant of areas" C = (q — c) A q is invariant 
along the solutions of (6.1). It is "projectivized" in C = c A ir 6 /\ 3 V, where 
7r = q Aq £ /\ 2 V is the projective impulsion: C is constant along the projective 
trajectories defined by the projective force field (6.2). On the other hand, the 
linear transformations of the group G introduced in Lemma 6.3 preserve C. 

Proof. Using the equation of dynamics (4.4), we see that C is invariant along a 
projective trajectory: d q C = c A dqir = c A / = ^(q)c Aq/\c = 0. Let jeG. We 
decompose it arbitrarily: tt = q A q' . As g(q) = q + 7c and g(q') = q' + 7'c for 
some (7, 7') £ E 2 , we have C~cAqAq' = cA g(q) A g(q'). 

Two-dimensional dynamics. The dynamics of (6.1) and of the projectivized 
versions (6.2) and (6.4) are essentially 2-dimensional. We will content ourselves 
with a study of the case n + 1 = dimU = 3. Hypothesis (6.3) corresponds 
to what is called in [A1S] a Jacobi attractor. Some conditions are given there 
ensuring that under this hypothesis, (6.1) possesses open sets filled by periodic 
orbits. Let us show how we arrive at the conclusions of [A1S] about the dynamics 
around a Jacobi attractor by an elementary study of the action of the group G. 

6.5. Lemma. The group G of Lemma 6.3 is abelian and canonically isomorphic 
to the additive group [c]° C V* , the annulator of the "vertical" line [c]. Assume 
dimU = 3. The domain 14 = V \ [c] is foliated by vertical half-planes with 
boundary the vertical line [c] that we simply call the leaves. Let us fix such a 
leaf L C V and call L C S(V) the corresponding object; L is invariant by the 
action of G. We call Lq C T~ X IA the set of projective states with position in 
L and given constant of areas C ^ 0; then G acts simply transitively on Lc, 
giving to Lc the structure of an affine plane, and to the submanifold of states in 
T~ X U with constant of areas C the structure of a principal G-bundle with base 
the one-dimensional manifold S{V/{c\). 
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Proof. The elements of G have the form Id + ui ® c £ V™ ® where w € V* 
satisfies (u>,c) = 0. The first statement follows easily. Let us fix a leaf L, and 
choose linear coordinates (x, y, z) of V such that c = (0, 0, 1) and L is y = 0, 
a; > 0. To give an element of Lc> we give a g with y(q) =0 that we normalize 
with the condition x(q) = 1, then we give a q such that g A q is the projective 
impulsion. We choose q with x(g) = 0, and we must take y(q) = (dxAdyAdz, C). 
After these normalizations only the z-components of q and q are free. But 
g 77 ' = Id + (jdx + j'dy) <8> c S G sends (g, g) onto (q + 7c, g + j'c). 

6.6. Proposition. Consider Problem (6.4), defined on the domain U = V \ [c]. 
Assume moreover dimF = 3, fix an integer k > and a constant of areas 
C ^ 0. Consider a "source" leaf L C V and "target" leaf if C V. Any initial 
condition in Lc (defined in 6.5) defines a projective trajectory that cuts the leaf 
K infinitely many times "afterwards" . The map from Lq to Kc that associate 
to the initial condition in Lc the state at the intersection with K after k turns 
commutes with the action of G. 

Proof. The domain U C S{V) is topologically the sphere minus two points, the 
"poles" , and a trajectory with C ^ cuts the projected leaves, the "meridians" , 
transversally. We prove first that any such trajectory is made of an infinity 
of "loops" around U. If it was not the case, the trajectory would stop at a 
pole, and by monotonicity of the "longitude", would do it tangentially to a 
meridian. But in a small sector near this meridian, the field of forces is close to 
the field of forces in a Kepler problem, for which such behavior is impossible. 
So the trajectory must turn indefinitely in the past and in the future. The 
commutation statement is Lemma 6.3. 

6.7. Eccentricity vector. For any orbit with C ^ 0, there is a unique element 
g G G that sends a state in Lq to the "next" intersection of the orbit with L. 
If g ^ Id then the space of orbits with same C is a cylinder, isomorphic to the 
quotient of G by the subgroup generated by g, and these orbits are not closed. If 
g = Id we are in what we called in [A1S] and [All] the Jacobi-Darboux case. The 
orbits are closed. In term of the G-principal bundle described in 6.5, Problem 
(6.4) defines a local trivialization of the bundle, and in the case g = Id this local 
trivialization has no monodromy. 

This means that we got a first integral for the states with given C. When g = Id 
it takes value in one of the fibers Lc ■ If we choose a reference orbit with this 
angular momentum 2 , the value of the first integral may be identified to the 
element of G that sends the reference orbit onto the present orbit. The element 
of G are covectors (see 6.5). The covectorial value of the first integral is related 
to the well-known eccentricity vector. 

We have all in mind that to a first integral should be associated a symmetry. 
Nothing excludes that in the same problem different "associations" coexist. 
Here we can speculate that to the action of the abelian symmetry group G is 
associated a "moment" that may be an element of the Lie algebra of G, and 

2 In the Kepler problem, for which g = Id, the natural choice is the circular orbit. 
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may be the eccentricity vector. 

If we pass to the particular case of the Kepler problem, which has a natural 
Hamiltonian structure, this association coexists with the Hamiltonian associa- 
tion, where the "moment" is the eccentricity vector multiplied by the square 
root of the given semimajor axis, and the Hamiltonian action is the rotation of 
the Bacry-Gyorgyi parameters (as defined in [Sou] or [Cor]). 

Up to now we gave to C a fixed non-zero value. Is the dynamics similar for 
the other non-zero values of C? The Jacobi problem (6.3) in affine dynamics 
possesses an obvious invariance by dilation (q — c,q) (A(q — c), A -1 / 2 */), which 
does not preserve the time parameterization. From the projective view point, 
this transformation is defined after the choice of an affine screen, and it does 
not commute with G. Nevertheless it extends to U and preserves the leaves. 
If we call (z, y, z) — (z(q), y(q), z(q)) the coordinates after the normalization in 
the proof of 6.5, the map is (z,y, z) ^ (\~ 1 z, A _1 / 2 i), A~ 3//2 i). This map may 
be used to establish that if the orbits are closed for a given non-zero value of C, 
they are closed for any non-zero value of C. 

7. The Kepler problem 

If (6.1) is (1.4) with [3 = —3, i.e. q = — \\q — c||~ 3 (<7 — c), the properties obtained 
in the previous section apply. The field of projective forces extends to V \ [c]. 
Its expression is: 

/(«) = IMI" 3 «Ac, (7.1) 

where |j.|| 2 is a positive quadratic form on V, satisfying ||c|| 2 = and having 
only [c] as a direction of degeneracy. 

7.1. Proposition. Let H be the group of linear transformations of V that 
preserve the degenerate quadratic form ||.|| 2 and fix c. Let dim V = n+ 1. The 
dimension of H is n(n + 1) /2. An element of H sends any projective trajectory 
solution of (7.1) onto another solution, and the constant of areas C is preserved. 

Proof. Let us choose a vectorial hyperplane W C V that does not contain c, and 
I E H. The image l(W) of W is a hyperplane that does not contain c. But the 
group G of Lemma 6.3 is a subgroup of H . We know that it possesses a unique 
element g that sends W on l(W). Now g^ 1 o I = k is an isomctry of W for the 
induced Euclidean form. Conversely any g G G composed with any isometry k 
of W is an clement of H. Then dimi? = dimO(VK) + dimG = n(n - l)/2 + n. 
The last statements come from the invariance of / under the action of H, and 
from the invariance of C stated in 6.4. 

Problem (7.1) may be seen as the "abstract" Kepler problem. The usual Kepler 
problem in the Euclidean space, and its 19th century generalizations to spaces of 
constant curvatures, appear as "materializations" obtained by choosing different 
screens. We will come back to these screens in a forthcoming work 3 . They all 
remove part of the symmetry of the abstract problem. 

3 See [BoM] for these generalizations. The remark that these different Kepler problems are 
related by central projection and change of time is due to Appell (see [Apl], p. 158). 
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But there is also a canonical screen: the screen \\q\\ — 1, the "unit cylinder" of 
the degenerate quadratic form. The group H is a symmetry group for this cylin- 
der. In a sense \\q\\ = 1 is more symmetric than a usual Euclidean cylinder, be- 
cause it admits the subgroup G of 6.3, whose elements are called "transvections" , 
as a group of symmetry. But in another sense it is less symmetric, because the 
translations in the c direction are not in the symmetry group. Restricting Prob- 
lem (7.1) to this screen we "materialize" the abstract Kepler problem without 
removing any symmetry. 




7.2. Figure. Construction of a Keplerian cone. 

As usual we think of c as the "vertical" direction. The field of forces induced 
on the screen by restriction of (7.1) is vertical and constant: the screen function 
is h(q) — \\q\\, so on h{q) = 1 = (dh,q) the restricted force fh = dh\f = 
(dh, q)c — (dh, c)q is c. The next proposition will describe the orbits. But we 
begin with the free motion, i.e. the case of a zero force. 

We know that a free motion describes the intersection of the cylinder = 1 
with a vectorial plane. This does not look like a free motion on a usual cylinder, 
but one should not forget that in screen dynamics, the reaction, that maintains 
the particle on the screen, is "radial" and not normal (see §4.6). The screen 
dynamics equation for the free motion is q — Xq, where A is imposed by the 
constraint \\q\\ = 1. If we project vertically this equation, i.e. if we project 
by the canonical projection V — ► V/[c], q <— * q c , the equation remains of the 
same form q c = Xq Cl which implies that the "horizontal motion" is uniform (the 
reaction being normal to the circle = ||g c || = 1). Let us call 9 the constant 
angular velocity of this motion. Then A = — 9 . 
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7.3. Proposition. The screen dynamics of the Kepler problem (7.1) restricted to 
the cylindrical screen ||g|| = 1 is as follows. If C = the particle describes the 
vertical lines (with direction c) in a uniformly accelerated motion. If C ^ the 
particle moves along the intersection of the cylinder with a non-vertical affinc 
plane, and does it in such a way that the canonical angular variable on the 
cylinder increases or decreases uniformly. 

Proof. The equation is q = Xq + c. The constant of areas is C = c A q A q. As 
in the case of a free motion, 9 is a constant proportional to C and A = — 9 . 
If we set q = 9~ 2 c + go, the equation becomes <jo = Xqo, equation of the free 
motion. The conclusion follows from the study of the free motion. The affinc 
plane passes through the vector 9 c. 

The intersection of the affine plane with the cylinder is a conic section. In a 
projective trajectory of the Kepler problem the ray-particle describes a quadratic 
semi-cone. The most famous result of Newton's Principia, that the orbits of the 
usual Kepler problem arc conic sections, is an obvious corollary. The above 
deduction of this famous result is indeed a kind of geometrization of the well- 
known deduction where one uses as variables the inverse of the distance from 
the particle to the center and the polar angle 9. 

7.4. About the parameterization of the paths. A projective trajectory possesses 
a field of tangent projective impulsions. Does it possess also, canonically, a 
field of tangent velocity vectors? In another way (see 3.3), is it canonically a 
usual trajectory? The answer is no for the uniform motion, but yes for the 
Keplerian motion: the cylindrical screen is canonical. It converts canonically 
the field of projective impulsions into a field of velocity vectors. But still the 
general "chronological" structure is the tangent field of projective impulsions, 
even if in particular situations some classes of parameterizations may appear as 
concomitant to the given field of projective forces. 

Let us visualize the "web" of all the projective lines drawn in the projective 
space. A point and a tangent direction uniquely determine a line. The "pro- 
jective geometry of paths" , a classical subject founded by the papers of Weyl 
in 1921 and Eisenhart in 1922, describes some possible "deformations" of this 
web, where a point and a tangent direction still determine a curve. It is possible 
to introduce in this theory the concept of projective impulsion. A projective 
impulsion is more than a tangent direction. However changing the length of the 
projective impulsion without changing its direction we do not change the curve. 

An extension to the projective geometry of paths has been investigated by 
Tabachnikov (see [OvT]): the "projective billiards". Here also the notion of 
projective impulsion could be introduced. We could add to the model a field of 
projective forces, thus "deforming" again the web of curves. Then the curves 
with same initial direction but distinct initial projective impulsion separate into 
distinct trajectories. 

There is a deep relation between the existence of interesting rcparameterizations 
and the intcgrability of a system. The discoveries by Matveev, Topalov [MaT] 
and Tabachnikov [Tab] on the geodesic flow on the ellipsoid in any dimension, 
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and their subsequent investigations on the subject are remarkable (see [Tal]). 

To come back to our initial question, let us mention that several works on the 
projective geometry of paths discuss parameterizations or classes of parametriza- 
tions of the paths. Works by T.Y. Thomas and by Whitehead are often men- 
tioned. However, as far as we know, these constructions are not canonical, and 
involve more structure than what is strictly necessary to define the paths. 

In forthcoming papers we wish to insist on the notion of projective impulsion. 
We can do this by presenting the focal theory of conic sections in pseudo- 
Ricmannian spaces of constant curvature, in such a way that their role in the 
Kepler problem is deduced very simply. We will also explain some results on the 
first integrals of the uniform motion, which are useful in the study of projective 
dynamics. Our projective impulsion point of view clarifies a classical literature 
on the subject culminating with a result of Nijenhuis 4 . 

8. Bibliographical Notice 

Halphen's remark. In 1878, Georges- Henri Halphen explained some of the ideas 
we presented in 1.1 and 1.3, beginning with these words: "On me permettra, 
en dernier lieu, d'appeler l'attcntion sur ce fait curieux..." This is, as far as we 
know, the first time a projective transformation is considered in dynamics. 

Halphen introduced a projective transformation in order to synthesize his answer 
to a question raised by Bertrand. Bertrand wanted to know what are the fields 
of forces in an affine space, depending only of the position, such that the particle 
describes a conic section, whatever be the initial condition. 

In 1877, a week after the publication of his question, Bertrand announced that 
he had just completed a first step toward the solution: he proved that the 
force had to be central. Darboux attended his lecture at the College de France 
and immediately completed the solution. Soon after Halphen presented his 
independent resolution 5 . 

Darboux and Halphen found two classes of solutions, and Halphen found a 
wonderful description of these classes: they are the field of Newtonian attraction 
from a center and all its affmc transforms and the field of linear forces (i.e. 

4 [Nij], For shorter proof see [Thl]. Many thanks to Bartolome Coll who introduced me to 
this literature. 

5 The references are the volume 84 (1877) of Comptes rendus de l'Academie des Sciences, 
pp. 673, 731, 760, 936, 939, and [Hal] for the remark on projective transformations. Actually, 
it is not clear that Bertrand's problem has ever been completely solved. An objection appears 
in a parenthesis of [Rou], §393 "unless the force at every point of that curve is infinite". The 
context of the original question makes clear that Bertrand looks for all the force fields in the 
plane. To deduce that the force field is central, we must use Bertrand's argument or Halphen's 
argument. Halphen's argument is 3-dimcnsional and does not work in the plane. Bertrand's 
argument does work and shows that the lines of force are straight lines. But these lines could 
be, for example, the tangent lines to a strictly convex domain. The force field would only be 
defined out of this domain. In [Imc] there is the claim of a complete solution of Bertrand's 
question, but the argument is clearly wrong. If we start at p. 32 with a force field (X, Y) 
which is central and satisfies Bertrand's requirement, but whose center is not the origin of the 
coordinates, we conclude at p. 36 that the origin is also a center of force, which is absurd. 
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q = —q) from a center and all its projective transforms. In fact we should accept 
imaginary parameters in the transformations to allow changes in the signature, 
or keep everything real and complicate a bit the statement. But our simplified 
statement raises immediately the question: why not the projective transforms 
of the Newtonian attraction also? And the similar but easier question about 
the linear force? The answer to the first question is: the Newtonian attraction 
is invariant by the projective transformations. This is the starting point of 
projective dynamics. 

Some earlier works by Newton and Hamilton describing the second class were 
mentioned by Glaishcr (see also [Tis] p. 42). But apparently they were not 
followed by a remark on the projective transformation. 

The articles [App] and [Apl] by Paul Appell on the central projection present 
several results, including some characterization of the projective transformations 
within a larger class of transformations acting on the space-time. Kasner later 
improved this characterization and extended it with De Cicco. 

Appell's work was mainly continued through a question suggested by Goursat 
(see the end of [App]) and the papers of Painleve, Rene Liouville and Levi- 
Civita dedicated to it (see also [Lut]). These works remained quite far from the 
original remark by Halphcn. Also in my sense the problem stated by Goursat 
does not respect the "non-quadratic" or affine character of equation (1.1). As 
far as I know, all the related and subsequent works involve a metric, even if 
T.Y. Thomas happily joined the question to the projective geometry of paths, 
and stated it using a connection (see [Tho] and [LiA]). 

Acknowledgements. We wish to thank Alain Chenciner, Mauricio Garay and 
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